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ABSTRACT 

EUV scatterometry is performed on 3D patterns on EUV lithography masks. Numerical simulations of the 
experimental setup are performed using a rigorous Maxwell solver. Mask geometry is determined by minimizing 
the difference between experimental results and numerical results for varied geometrical input parameters for 
the simulations. 

Keywords: EUV scatterometry, optical metrology, 3D rigorous electromagnetic field simulations, computational 
lithography, finite-element methods 

1. INTRODUCTION 

Extreme ultraviolet (EUV) lithography at a wavelength of about 13 nm is expected to replace deep ultraviolet 
(DUV) lithography for fabrication of integrated circuits with minimum feature sizes (critical dimension, CD) as 
small as 22 nm or below. With decreasing feature sizes tolerance budgets of mask pattern dimensions get tighter. 
Metrology of such structures has to be performed for characterization and process control. 

In previous works it has been shown that EUV scatterometry is a fast and robust method for characterizing 
masks with ID-periodic patterns (line masks). I^^Sl In this contribution we extend our research to the characteri- 
zation of 2D-periodic patterns (array of contact holes) . 

In our method we compare the distribution of EUV light scattered off the sample to results from numerical 
simulations of the setup with different sets of geometrical parameters. Mask geometry is then determined by 
optimizing geometrical input parameters of the numerical simulations such that differences between measurement 
and simulation are minimized. 

This paper is structured as follows: The experimental setup of the EUV reflectometer at PTB is presented 
in Section [2] The used model of the investigated mask and the numerical method are described in Section [31 
Section H] shows results on 3D geometry reconstruction. 
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Figure 1. Left: Scheme of scatterometry measurements: The detector angle /3 is scanned at a fixed incidence angle a. 
Right: Photograph of the experimental setup. 



2. EXPERIMENTAL EUV SCATTEROMETRY RESULTS 

The data presented here were measured at the EUV reflectometry facility of PTB^"^ in its laboratory at the 
storage ring BESSY II. The PTE's soft X-ray radiometry beamlin^ uses a plane grating monochromator which 
covers the spectral range from 0.7 nm to 35 nm and which was particularly designed to provide highly collimated 
radiation. For this purpose it uses a long focal length of 8 m in the monochromator, and the focusing in the non- 
dispersive direction is provided by a collecting pre- mirror with a focal length of 17 m. We achieve a collimation 
of the radiation in the experimental station to better than 200 /.trad and the scatter halo of the beam can be 
suppressed to below 10~^ relative intensity at an angle of only 1.7mrad to the central beam. 

The measurement scheme for angular resolved scatterometry of periodic structures is presented in Figure [1] 
{left). Usually, the sample is set at a fixed angle to the incoming radiation and the detector is moved. Figure [1] 
(right) shows a photograph of the experimental setup: The probed photomask is mounted at the center of the 
(red marked) x-y-z-coordinate system. It is rotated such that the periodicity vectors of the periodic patterns on 
the mask are in x- and y-directions. The EUV probe beam is incident in the y-z-plane. The detector can be 
freely moved to capture reflected x- and y-diffraction orders. 

The presented experimental data are obtained from periodically structured areas on an EUV photomask. 
The structure for the measurements was a 1.4 mm by 1.4 mm large field with a 2D structure of quadratic contact 
holes with nominal 300 nm CD at 600 nm pitch in x- and y-direction. The measurement method used was angular 
resolved scatterometry. The scan range of the detector angle /? was from 1° to 25°, see Figure [1] The —8*'' 
to -1-6*'' diffraction orders were measured in y-direction with respect to the sample using the movement of the 
detector in the plane of specular refiection. For the measurement of diffraction in the out-of-plane direction, 
the perpendicular movement of the detector was used.- Thus we also measured from the —2""^ to -1-8*'' order in 
x-direction. As the EUV radiation is resonant with the underlying multilayer refiective coating, the measured 
diffraction intensities strongly depend on the actual wavelength. We measured at centre wavelength of the 
resonance peak of 13.54 nm, see Figure [31 Figure [5] shows a scatterogram for a measurement of the contact hole 
array in two different graphical representations. 
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Figure 2. Scatterogram of an EUV mask with a contact hole array. Left: Measured diffraction intensities (logioil]^^ )) 
as function of diffraction index in x- and t/-direction, A^x, Ny in a pseudo-color representation. Right: Same data as a 2D 
plot on a logarithmic scale. The dashed line indicates the threshold of 8 ■ 10^^ used in the evaluation. 
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Figure 3. Reflection from the blank multilayer mirror (without buffer, absorber and ARC): Experimental results (o) 
numerical results (x) for multilayer parameters as defined in Table [T] 
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3. MASK MODEL AND SIMULATION METHOD 
3.1 Multilayer model 

For characterizing the multilayer stack we have performed measurements of reflectivity off the blank multilayer 
mirror as a function of incident wavelength in a spectral range from 13 nm to 14.3 nm. The multilayer stack 
consists of 40 A/2-multi- layers of molybdenum (Mo) and silicon (Si). It further contains a silicon capping layer, 
finished with a thin oxide layer, and it is placed on a Si02 substrate. To model the A/2-multi-layers we have 
chosen a four- layer stack consisting of a Mo-layer, a MoSi2-layer (interlayer), a Si-layer, and a second MoSi2-layer. 
A schematic is shown in Figure |4] ( left) . We have fitted the layer thicknesses of the blank mirror using the IMD 
multilayer programme,'^ using also the therein provided material definitions. The resulting reflectivity spectrum 
is display in Figure |3] together with measured data. The obtained layer thicknesses are shown in Table [TJ These 
values are used in the subsequent 3D simulations. For the thicknesses of the absorber stack (buffer, absorber 
and anti-reflection, ARC, layers) we have chosen the nominal values of the EUV mask (see Table [T]). 



material 


height 


n 


k 




vacuum 


inf 


1 







ARC 


^ARC = 12 nm 


0.9474 


0.0316 


Rz — 5 nm 


absorber 


ft-a — 55 nm 


0.9255 


0.0439 




buffer 


/lb = 10 nm 


0.9735 


0.0131 




oxide 


/lox = 1-35 nm 


0.9735 


0.0131 




capping 


^cap = 8.19 nm 


1.0097 


0.0013 




multi-layer (interlayer) 


hi = 1.48 nm 


0.9681 


0.0044 


40 layers 


multi-layer (Mo) 


/lino = 2.99 nm 


0.9206 


0.0065 


40 layers 


multi-layer (interlayer) 


hi = 1.48 nm 


0.9681 


0.0044 


40 layers 


multi-layer (Si) 


hsi ~ 1.18 nm 


1.0097 


0.0013 


40 layers 


substrate 


inf 


0.9735 


0.0131 




lateral dimensions 


Ci^bottom,x 


300 nm 






(nominal) 


C'-C'bottom.y 


300 nm 










600 nm 








Py 


600 nm 








-Rxy 


Onm 








sidewall angle a 


90 deg 






illumination 


angle of incidence 


6deg 








wavelength Aq 


13.54 nm 






polarization 




r,0,0) 





Table 1. Parameter setting for the EUV mask simulations (compare Fig. [4]): Mask geometry parameters (layer heights h^, 
sidewall angle a, corner rounding radius R^), material parameters (real and imaginary parts of the refractive index, n and 
k, at a wavelength of 13.54 nm), illumination parameters (in-plane angle of incidence, vacuum wavelength, plane- wave 
electric field polarization). 



3.2 Lateral layout: array of contact holes 

The investigated mask holds several different patterns. Here we concentrate on analysis of a single of these 
patterns: an array of contact holes. Figure |3] (right) shows a schematic of the lateral layout. A pattern with a 
pitch of Px — Py — 600 nm is investigated. Nominal values for the (bottom) CD's are CD^ = CDy ~ 300 nm. 
We assume a sidewall angle of the absorber stack, a, and a constant (vertical) corner rounding of the stack, R^. 
Free parameters in our investigation are the contact hole critical dimensions CD^, C'Dy, the sidewall angle a, 
and the lateral corner rounding radius, i?xy 

3.3 Simulation setup 

For rigorous simulations of the scattered EUV light field we use the finite-element (FEM) Maxwell solver JCM- 
suite. This solver incorporates higher-order edge-elements, self-adaptive meshing, and fast solution algorithms 
for solving time-harmonic Maxwell's equations. Previously the solver has been used in scatterometric investiga- 
tions of EUV line masks (ID-periodic patterns). ^^1^ Recently we have reported on rigorous electromagnetic field 
simulations of 2D-periodic arrays of absorber structures on EUV masks.l^ This report contained a convergence 
study which demonstrates that highly accurate, rigorous results can be attained even for the relatively large 3D 
computational domains which are typically present in 3D EUV setups. 

Briefly, the simulations are performed as follows: a scripting language (Matlab) automatically iterates the 
input parameter sets (CD^, CDy, a, Rxy)- For each set, a prismatoidal 3D mesh is created automatically by 
the built-in mesh generator. Then, the solver is started, postprocessing is performed to extract the diffraction 
order efficiencies, and results are evaluated and saved. Figure [5] shows a graphical representation of a 3D mesh. 
As numerical settings for the solver, finite elements of fourth-order polynomial degree, automatic settings for the 
transparent boundary conditions, and a rigorous domain-decomposition method (to separate the computation 
of the light field in the multi-stack mirror from the computation in the structured absorber region) are chosen. 
Together with a relatively coarse lateral spatial discretization and a fine spatial discretization in ^-direction this 





Figure 5. Mesh of the structured absorber geometry. The discretizations of the adjacent vacuum region and of the 
multi-layer stack are not displayed. 

setting yields discrete problems with between one and two millions of unknowns. These problems are solved by 
direct LU factorization on a computer with extended RAM memory (about 40 GB) and several multi-core CPU's 
(total of 16 cores), with typical computation times of about 30 minutes (per parameter set). From comparison 
to a detailed convergence study on 3D EUV simulation^iSI estimate that we achieve relative accuracies better 
than 1 % for all diffraction orders with intensities greater than 10~^ (relative to the intensity of the incoming 
beam). 

4. RECONSTRUCTION OF GEOMETRY PARAMETERS 
4.1 Simulation results 

The simulated scatterogram of the EUV mask pattern, where all parameters are set to nominal values, is shown 
in Figure H) The left part of the Figure shows the diffraction intensities on a logarithmic, pseudo-color scale. 
The right part of the Figure displays the same simulation results in a 2D graph, together with results from the 
measurement. Only results with measured relative intensities larger than 8 • 10~^ are displayed (c/., dashed line 
in Figure [21). 

Each simulated diffraction order intensity, I^^n^ , is attributed a relative deviation: 




iVx diffraction orders 



Figure 6. Scatterogram of a contact hole array on an EUV mask. Left: Simulated diffraction intensities {logio{I%^jsf^)) 
as function of diffraction index in x- and y-direction, A^x, Ny in a pseudo-color representation. Parameter setting: 
CDx = CDy — 300 nm, a — 90deg, i?xy ~ Onm. The red, dashed line indicates the range of diffraction orders obtained 
experimentally. Right: Simulations (x, same data as left) and corresponding experimental values (o), intensities > fO~^ 
displayed on a logarithmic scale. The 2D array of orders is plotted row by row in y for subsequent orders in x. 



where /^^^ is the measured value and and A^y denote the index of the diffraction order in x-, resp., y- 
direction. 

A cost function, is attributed to the simulation parameters by summing over all diffraction orders iVx, Ny, 
with I'^^pf > 8 * 10~^/o, where Jo is the intensity of the illuminating plane wave: 

5P'^{CD^, CDy, a, R^y) - {^InT,n,{CD^, CDy, a, R^y) /n„,y . 

Ny 

Here, rim denotes the number of diffraction orders taken into account, and the zero**^ diffraction order is omitted. 
4.2 Reconstruction results 

We have performed a series of simulations with varied parameter sets {CD^, CDy, a, Rxy)i recorded the cor- 
responding diffraction patterns and computed the corresponding cost functions / deviations W™. Figure [7] 
{left, center) shows how the cost function varies with CD^ and i?xy, in the investigated parameter ranges. We 
find the lowest cost function, i.e., the best correlation between experiment and simulation for the parameter set 
CDx = 310 nm, CDy — 307 nm, a — 88deg, R^y — 90 nm. The measured and simulated diffraction spectra 
for this parameter set are displayed in Figure [7] (right). The agreement of most diffraction orders intensities is 
very good. We expect that the parameter values obtained with this scatterometric method are close to the real 
parameter values. 

However, the fact that the differences between measurement and simulation are larger than both, expected 
experimental measurement uncertainty and numerical errors, suggests the following advancements: 

(i) More free parameters to be considered in the model (e.g., layer thicknesses /iarc, ^a, ^b)- 

(ii) Larger parameter ranges to be scanned to avoid possibilities of optimizing to a local minimum. 

(iii) A more refined geometrical model to take into account effects of edge roughness j^^i^l surface roughness, 
corner rounding at the bottom corners, etc. 

(iv) Reconstruction results to be validated by comparison to results obtained with other measurement methods, 
like scanning electron microscopy (SEM) or atomic force microscopy (AFM). 

Currently, computation times for the 3D simulations on large computational domains limit the applicability 
when many free parameters are taken into account (unless large computational power is available). For real-time 
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Figure 7. Reconstruction results. Le/t; Dependence of cost function on CDx- Center: Dependence of cost function on 
i?xy Right: Simulated diffraction intensities (x) for parameter setting with highest correlation to experimental data (o). 



reconstruction of using the same methods we plan to apply the reduced-basis method (RBM) which allows to 
decrease computation times for rigorous simulations of parameterized simulation setups by several orders of 
mag nitude.lSlIll 

5. CONCLUSION 

Scatterometric measurements of 2D periodic patterns on EUV masks have been performed using the X-ray 
radiometry beamline of an electron storage ring. Rigorous simulations of the measurements have been performed 
using a finite-element based Maxwell solver. Very good agreement between experimental results and simulation 
results has been achieved. Parameter reconstruction has been demonstrated. 

For the future we plan to perform measurements of further and more complex patterns, comparison to AFM 
and SEM results, and application of reduced basis simulation methods. 
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